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Trapped ions arranged in Coulomb crystals provide us with the elements to study the physics of 
a single spin coupled to a boson bath. In this work we show that optical forces allow us to realize 
a variety of spin-boson models, depending on the crystal geometry and the laser configuration. We 
study in detail the Ohmic case, which can be implemented by illuminating a single ion with a 
travelling wave. The mesoscopic character of the phonon bath in trapped ions induces new effects 
like the appearance of quantum revivals in the spin evolution. 

PACS numbers: PACS 



The problem posed by a two-level system interact- 
ing with a bath of harmonic oscillators, known as the 
spin-boson model, appears in condensed matter, atomic 
physics, and quantum information processing. It is of 
fundamental importance, since it represents a paradigm 
for the study of quantum dissipation and the quantum- 
to-classical transition [H, 0, The spin-boson model 
displays nonperturbative features of many-body physics, 
like the localization transition at a critical dissipation 
strength. For the special case of an Ohmic bath spec- 
trum, it may be mapped onto the celebrated Kondo 
Hamiltonian. Despite its fundamental importance, ex- 
perimental investigations into anything but the weak- 
coupling regime of the spin-boson model are still scarce. 
The localization transition has been observed in the re- 
lated Josephson junction systems Q, while typical solide- 
state two-level systems feature a coupling strength much 
below the critical treshold [5|. 

Trapped ions provide a clean experimental model sys- 
tem which is ideally suited for the quantum simulation 
of condensed matter problems @, 0| • As we show in this 
work, the spin-boson model appears here in a natural 
way, with a wide range of tunable parameters. The role 
of the two-level system is played by two internal levels of 
a single ion that is part of a Coulomb crystal. Since ions 
interact through the Coulomb repulsion, their motion is 
described by collective modes with a dispersion relation 
which depends both on the direction of vibration and on 
the trapping conditions. The internal level of the sin- 
gle ion can be coupled to the vibrational bath by means 
of optical forces [8(, in such a way that several differ- 
ent spin-boson couplings may be implemented. Trapped 
ions offer us a wide range of possibilities for observing 
the phenomenology of the spin-boson model, with the 
advantage that we can use well-developed experimental 
techniques for the preparation of both the initial spin and 
phonon bath states [9j. Besides that, the finite number 
of vibrational modes forms a mesoscopic environment. 
This allows us to study memory effects due to the finite 
size of the system, such as quantum revivals in the spin 



evolution, which are absent in the customary limit of a 
macroscopic bath. 

In this work we analyze the spin-boson models which 
appear in systems of trapped ions, and we show the fol- 
lowing results: (i) ID and 2D Coulomb crystals provide 
us with a variety of phonon spectral densities, ranging 
from sub- to super-Ohmic, depending on the ion crystal 
dimension and the laser configuration, (ii) The particular 
case in which the internal level of a single ion is coupled 
to a Coulomb chain by a travelling wave corresponds to 
the Ohmic spin-boson model. If the number of ions is 
large enough, this setup allows us to tune the interac- 
tion strength and recover the standard phenomenology 
of this model, such as the quantum phase transition to- 
wards localization at strong coupling, (iii) For time scales 
larger than a given revival time, finite size effects induce 
the reexcitation of the spin after an initial period of re- 
laxation (quantum revival), which can be observed in a 
wide range of parameters, including high temperatures 
of the phonon bath. 




FIG. 1: (color online). Scheme of the implementation of spin- 
boson models in a Coulomb chain. 

Experimental realization. — The standard spin-boson 
Hamiltonian is defined by 

H = l^ + yO'x+X] ( W "4 a n + y A n( a n +«")) ■ (!) 
n 

a a are the Pauli matrices, and lu u and a n are the en- 
ergies and annihilation operators of the bath phonons, 
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respectively (we set h = 1). Two hyperfine levels | — 1) 
and |1), separated by the internal energy u>o, correspond 
to the eigenstates of a z . We propose to experimentally 
realize the Hamiltonian |T]) by addressing a single ion in 
a Coulomb crystal in one of two laser configurations (see 
Fig. [T|): (i) State dependent dipole force. In this setup, 
a laser induces the tunneling term in Eq. ([1]), in such a 
way that the tunneling amplitude A and the bias e are 
set by the laser's Rabi frequency, fit, and the detuning, 
u)q — wl, respectively. The spin-phonon coupling is in- 
duced by subjecting the selected ion to an off-resonant 
standing wave, which creates the state dependent dipole 
potential V(z) — Vq cos 2 (kz + <fi) a z . k and <j> are the 
wavevector and phase of the standing wave lasers, re- 
spectively. The operator z is the ion's coordinate in the 
direction of the optical force relative to its equilibrium 
position. It is readily expressed in terms of phonon op- 
erators, z = J2 n MnZn(a n +<4), where z„ = l/y/2muj n , 
with m the ion mass, and M n is the amplitude of each 
vibrational mode n at the ion. In the Lamb-Dicke regime 
(kz <C 1) we can expand V(z), 



H sw = \a z (V(0) + Fz + Qz 2 



(2) 



V(0) is just a shift of the bias e, and by choosing 
appropiately one can retain the linear or the quadratic 
terms in (|5J) only. We focus on the case of a linear 
coupling (<fi — 7r/4, Q = 0), which leads to (Q} with 
X™ = FM n z n . (ii) Polaron coupling induced by a trav- 
elling wave. In this setup a single laser creates both the 
spin-boson coupling and the tunneling term in Eq. (TT|). 
One of the trapped ions in the chain interacts with a trav- 
elling wave, such that the coupling is given, in a frame 
rotating with wl, by 
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Hamiltonian ([3]) can be recast to fit the expression (JTJ) by 
applying the canonical transformation U = e -( l / 2 ) kzry ^ 
which yields A'™ = kz n A4 n uj ni A = Ql, and e = u>q — u>l. 

Spectral density. — The properties of the spin-boson 
model are determined by the spectral density, J(uj) — 
,r Er ( =i^n — w„). The model shows a rich phe- 
nomenology that has mostly been analyzed for the case 
of an infinite gapless bath, such that J(u>) a w s . In 
this work we focus on describing experimental conditions 
that lead to this situation. Two main cases are rele- 
vant, namely, ID and 2D Coulomb crystals. Although 
finite size effects will play an important role, and will 
be discussed below, we will first revisit the thermody- 
namical limit, N — ► oo, to understand some of the main 
features. The axial (ID) or in-plane (2D) modes can 
be considered approximately gapless, since the minimum 
energy is the global trapping frequency along the crys- 
tal, which has to be small enough to guarantee the sta- 
bility of the latter. For a ID Coulomb chain of ions 



with equally spaced equilibrium positions, the low energy 
spectrum of axial vibrations TO] is approximately given 

log(^)) 1/2 



by = (1 - f log (^)r' , where v = y/30/2, 

lu z is the axial trapping frequency, and f3 = 2e 2 /(mw 2 do), 
with do the mean distance. In order to get a qualitative 
description, we will retain the linear term only (see Fig. 
[2] (a)), and get the following spectral densities for each 
coupling: 



J™(w) - {l/v)M 2 {Fz Q fu-\ 

J™(U) - (1/V)M 2 V 2 LJ = 27TC* u, 



(4) 
(5) 



where M is the mode function of axial phonons at the ad- 
dressed ion, r) — fc/y / 2mtJ7 is the Lamb-Dicke parameter, 
and we have reexpressed J^(uj) as a function of the di- 
mensionless dissipation strength a [lj . Finite size effects 
will modify the spectral density of a Coulomb crystal. 
Nevertheless, Eqs. (j4} and (O allow us to determine the 
character of the phonon bath, as well as the scaling of its 
properties with experimental parameters. This is shown 
by the comparison with exact numerical calculations in 
Fig. O Finally, in a 2D Coulomb crystal, ions arrange 
themselves in a triangular lattice ll|. There, the low- 



est energy vibrational modes also show a linear disper- 
sion relation |12[ , such that the corresponding spin-boson 
models have an algebraic spectral density with exponents 
s = [l3| and s = 2, in the cases of interaction with a 
standing wave, or a travelling wave, respectively. 
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FIG. 2: Phonon bath properties of a Coulomb chain with 
iV = 50 ions, (a) Dots: Axial vibrational spectrum. Contin- 
uous line: Approximation given by retaining the linear terms 
in lj^ d . (b) Spectral density corresponding to spin-boson cou- 
plings induced by a standing or a travelling wave. The contin- 
uous lines are obtained with the exact vibrational spectrum, 
by substituting the delta function in J(ui) by a constant func- 
tion in each interval (w n ,w n +i), of height (uj n +i —cjn)" 1 . The 
dashed lines are the approximations of Eqs. (2} and (|5)l . 

Trapped ion Ohmic model. — In the following we ana- 
lyze in detail the spin-boson model induced by a travel- 
ling wave interacting with the central ion of a Coulomb 
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chain, which leads to ([5]). Since this case corresponds to 
Ohmic dissipation, it shows a rich phase diagram as a 
function of the temperature, T, and a Assume that 
at t < the coupling ([5]) is off, the phonon bath is in 
thermal equilibrium, and the internal state is prepared 
in |1) by using standard optical pumping techniques. We 
focus on the evolution of the system at time t > 0, upon 
suddenly switching on the laser. In particular we calcu- 
late P(t) = (er 2 ), to determine the evolution of the spin 
under the effect of the tunneling A [1]. The theoreti- 
cal description of this system can be addressed within 
the Non-Interacting Blip Approximation (NIBA), which 
consists of a non-Markovian evolution equation for P(t): 



P(t) = / K(t-t')P(t')dt', 
Jo 

K(t) = -A 2 Re(e lkz{T) e- lkz(0) ) 



(6) 



where the average is evaluated assuming a thermal state 
in the bosonic bath. This expression can be obtained by 
neglecting spin-bath correlations in a perturbative ex- 
pansion up to second order in the tunneling amplitude A 

HI- 

It is well established that it describes correctly the 
two limits of weak and strong dissipation, as well as the 
high temperature limit. 

The finite-size properties of the spin-boson model de- 
fined by Eq. ^ can be qualitatively understood by con- 
sidering the approximation that vibrational energies are 
equally spaced by a given energy Sto. Under this assump- 
tion the Kernel K(t) defined by ((6]) is periodic, and can 
be written as K(t) = ^ K(t — r„), r„ = nr rcv , where 
T rov = 2tt/Suj is the vibrational bath revival time. K(t) 
becomes equal to the Kernel of the continuous Ohmic 
model in the limit T rov — » 0. At short times it <C T rov ) 
the spin evolution is governed by K(t), and it shows a 
similar behavior as in the case of the continuum spin- 
boson model. The periodic structure of K(t) manifests 
itself at long times (r > t 1cv ) in the form of quantum 
revivals in P(t). The spin-boson model of trapped ions 
satisfies the condition of equally-spaced vibrational en- 
ergies only approximately. However one can define an 
average energy spacing Slo, which in the limit of large 
N can be estimated by 8lj vuj z Tt/N. Thus, in a real 
Coulomb chain we expect that r rev ~ 2-7T /8u> defines the 
time scale which separates the short-time regime in the 
evolution of P(t) coinciding with the evolution for a con- 
tinuous bath, and the long-time regime displaying the 
quantum revivals. 

Numerical solution. — To verify this assumption, we ob- 
tain the numerical solution of P(t) within the NIBA, us- 
ing the exact vibrational modes of a real ion chain. In 
order to compare the results obtained in the case of a real 
Coulomb chain with those predicted by the ideal Ohmic 
spin-boson model, we fit the exact low-energy spectral 
density to the form given by ([5]), and from there we ex- 
tract the dissipation strength a which describes the ion 



vibrational bath. The NIBA can be justified for a discrete 
phonon bath as long as the decay time of K(t) is much 
smaller than r rcv . At short times t <C T rcv , the validity of 
the NIBA in the continuous case implies that the aver- 
age in Eq. (JSJ) can be calculated by factorizing the total 
density matrix into the spin and phonon reduced density 
matrices. Due to the periodic structure of K(t), brought 
about by the discreteness of the vibrational bath, we con- 
clude that P(t) is related to the state of the system at 
times t — 7"iT rov - In the first revival (t T rov ), P{t) de- 
pends on the state of the system at short-times, where the 
spin-bath decoupling scheme, which leads to the NIBA, 
works. Thus, the first revival is well described within in 
the NIBA, and the argument can be easily extended to 
later revivals. We checked the NIBA for the finite-size 
bath vs. the Weisskopf-Wigner approximation at weak 
coupling, finding very good agreement. The numerical 
solution of the NIBA integro-differential equation is per- 
formed by means of the numerical method presented in 
[3] . We now discuss in detail the results for two different 
regimes, depending on the phonon temperature. 




FIG. 3: (a) Transition from the underdamped (thick line, 
A = Wloz, a = 2 ■ 10 -3 ) to the overdamped regime (thin line, 
A = 3uj z , a = 4 • 10 -3 ), in the high-temperature limit of a 
spin-boson model with N = 50 ions, T = 250w z . (b) NIBA 
time kernels corresponding to these two regimes. 

High temperature regime. — For sufficiently high tem- 
peratures, the kernel K(r) decays exponentially with a 
memory time r m = l/(2iraT) (we set Ub = 1)- In the 
continuum limit one finds two main regimes: r m A > 1 
(coherent oscillations) where Pit) oscillates with A and 
decays in a time r m ; and r m A < 1 (overdamped relax- 
ation) where Pit) decays with a rate T = A 2 r m . In the 
case of a finite Coulomb chain, K{t) shows an exponen- 
tial decay at short times, and additionally, an approxi- 
mate periodic structure at time scales r rev w 2tt/Suj (see 
Fig. [3] (b)). Fig. [3] (a) shows that the behavior of P(t) 
at short times clearly reveals the transition between the 
overdamped and underdamped regimes, as well as the 
quantum revivals at r rev . The revival effect can be un- 
derstood in terms of the perturbation of the Coulomb 
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FIG. 4: Evolution of the spin population in the Ohmic 
trapped ion spin boson-model, in a Coulomb chain with 
N = 50 ions, T = 0, A = lCkj z , and different dissipation 
strengths. 



chain created during the initial spin relaxation, which 
propagates along the chain, is reflected at the bound- 
aries, and returns to the selected ion, thus inducing its 
reexcitation. Interesting geometrical effects on the re- 
vivals may also be expected in 2D setups. Revivals in 
the high-temperature regime could be easily observed in 
experiments with trapped ions, since they require neither 
cooling to very low temperatures, nor high values of a. 

Low temperature regime. — For the Ohmic model in the 
continuum and scaling limits, the evolution of P(t) at 
T = is determined by the value of a, in such a way 
that there are three regimes to be considered [l|: a < 1/2 
(coherent oscillations), 1/2 < a < 1 (overdamped relax- 
ation), and a > 1, in which case dissipation impedes the 
decay of P(t) and the system becomes localized in the 
initial value of the effective spin [17| . Since the NIBA is 
known to reproduce the transition between these three 
regimes [Hj], we can use it to investigate whether the 
Ohmic spin-boson model in trapped ions shows the same 
transition. Our results are plotted in Fig. |4j where we 
show that the relaxation of P(t) shows the same qualita- 
tive features as in the standard Ohmic model, with the 
appearance of quantum revivals at long times. The ad- 
ditional localization of the spin state is clearly evident at 
values of a > 1, although a residual relaxation process 
still persists as a consequence of the discreteness of the 
bath. To quantify in more detail the transition to spin 
localization, we have calculated the initial decay rates at 
short times, as a function of a. Our results in Fig. Ela) 
show the slowing down of the spin relaxation with in- 
creasing a, as well as the effect of finite temperatures. 
Note that Fig. [5fa) corresponds to the finite-size version 
of the quantum phase transition to localization which is 
found in the thcrmodynamical limit of the phonon bath 
(where T would vanish above a = 1). 

The strong dissipation regime of the mesoscopic Ohmic 



FIG. 5: (a) Evolution of the decay rate in the strong dis- 
sipation regime. N = 50, A = 10o; z . (b) Axial trapping 
frequency of the trap that is required as a function of dissi- 
pation strength in units of the recoil energy, with different 
values of N. Eb, = k 2 /2m, with k the wavevector in @. For 
example, Be + , 2n/k = 300 nm, yields E R = 245 kHz. 

spin-boson model requires high values of a, and thus of 
rj 2 cx 1 /u> z . In the case of a chain with N = 50 ions, the 
transition to localization can be observed with u z of a 
few kHz. The axial trapping frequency that is required 
to realize a model with a given a decreases with N, see 
Fig. [5fb). This condition is difficult to meet in an exper- 
iment, due to the need to cool the Coulomb chain to low 
temperatures. However, it must be noted that ground 
state cooling is not required. Our calculations show that 
up to temperatures of the order of the axial trapping fre- 
quency the transition to localization can still be observed, 
whereas at higher temperatures it is smeared out. 

Outlook. — Let us comment on the possibility of imple- 
menting various interesting experimental situations other 
than the ones dicussed above. In particular, the cou- 
pling ^ allows us to implement a bath with l//-noise, 
a model which is relevant to the description of decoher- 
ence of solid-state qubits (both normal and supercon- 
ducting) Besides that, according to our discussion 
for the case of an off-resonant standing wave addressing 
the selected ion, it is possible to tune the spin-boson cou- 
pling in such a way that we implement spin-boson models 
with couplings quadratic in the bath coordinate. Again, 
these are of current concern in the description of su- 
perconducting solid-state qubits operating at the 'sweet 
spot' [13], and could be studied in much more detail in an 
ion chain model. Finally, the simultaneous coupling of 
several spins to the vibrational bath, by addressing sev- 
eral ions with lasers, would represent an implementation 
of a 'many spin-boson' model, where the interplay be- 
tween the phonon bath-mediated spin-spin coupling and 
dissipation could be studied. 
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